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Abstract 

Quantities like the fine structure constant at the pole of the Z boson and the 
anomalous magnetic moment of the muon are of profound importance for testing 
the Standard Model of elementary particle physics. Because these quantities are 
known with very high precision, deviations between experimental measurements and 
theory predictions for these quantities open a window for so-called "new physics", 
i.e. physics beyond the Standard Model. In this seminar new calculation techniques 
for the most unreliable part of the theory predictions, the hadronic contribution, are 
formulated, discussed and used for a prediction of these quantities. 
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1 Parameters for a precision test 



The running fine structure constant or QED coupling a at the scale of the Z° mass and 
the anomalous magnetic moment a M of the muon |I|, [^, |[] are two parameters which have 
gained great deal of interest during the last years. The reason for this interest is that these 
parameters can be determined experimentally with high precision. Theoretical predictions, 
therefore, can serve as precision tests of the Standard Model (SM) of elementary particle 
physics because deviations to the measurement opens windows for the appearence of so- 
called "New Physics". 

The QED coupling 

For the QED coupling, an accurate knowledge of a(Mz) is instrumental in narrowing 
down the mass window for the last particle of the Standard Model, the Higgs particle, 
which might have been already detected at the LEP collider at CERN at 115 GeV but 
which needs to be verified. It is worth to consider this parameter at the pole of the Z° 
because a huge amount of data have been taken there during a rather long period at the 
LEP I run. One could expect a rather precise experimental measurement at this point. 

The anomalous magnetic moment 

The above mentioned window seems to open already for the anomalous magnetic moment of 
the muon. In February 2001, the Brookhaven National Laboratory (BNL) reported about 
a precision measurement of the anomalous magnetic moment of the positive muon [EJ, 

c^ xp = (116 592 023 ± 151) x 10~ n (1) 

which had to be contrasted with the theory prediction. At the first moment, a deviation of 
roughly 2.6 standard deviations gave rise to suggestions for new-physics effects published 
during the last few months, including concepts of supersymmetry, leptoquarks, lepton 
number violating models, technicolor models, string theory concepts, extra dimensions 
and so on (cf. for instance Ref. |J). However, the situation changed suddenly when an 
error in the sign was found in the calculation of the light-by-light contribution, obtained 
by Tom Kinoshita and collaborators (see e.g. Refs. ||). This sign was later on corrected 
by the authors themselves M. The discrepancy was reduced again. New data sets taken 
in Novosibirsk again churn the water, but the physics community is split into two parts at 
least since these events. 

The hadronic contribution 

From the theoretical point of view the uncertainty in the determination of aiMz) and a M 
is dominated by the uncertainty of the hadronic contribution a had and aj^ ad . This talk, 
therefore, will deal with these contributions only. 
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2 The integral representation 



In this section I want to explain where the hadronic contribution comes from. Actually, 
the main ingredient to this contribution is the hadronic two-point correlator of the photon 
current, and the details and reasons of their origin are detailed in the following subsections. 



Hadronic contribution to the QED coupling 

The corrections Aa(q 2 ) to the QED coupling constant are given by the corrections to 
the photon propagator due to the calculation of a chain of inserted vacuum-polarization 
terms 
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— e 2 g 2 n(— q 2 ) = — A7raq 2 U(— q 2 ) and n(— q 2 ) is the two-point correlator 
ical function except for a cut along the real axis for q 2 < —4m 2 , 
oo. In order to understand what is meant by this, I take a simple 

Am 2 — q 2 . (3) 



where II 7 (— q 2 ) 

This function is an analytical function except for a cut along the real axis for q 2 < —Am 
and it vanishes for \q 2 
example, namely 
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This function takes an imaginary value 
for q 2 > Am 2 . The sign of the imagi- 
nary value depends on whether we ap- 
proach the real axis from the upper or 
lower half plane, as it is shown in the 
figure on the right hand side. This 
should be indicated by q 2 
s > Am 2 . So we obtain 
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The spectral density is finally given by 
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HM(g 2 ) = — 
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After having investigated the cut by introducing the dis- 
continuity, we can use Cauchy's theorem in order to ex- 
press the correlator function by the corresponding spec- 
tral density. For this we take the circle about the specific 
point z = q 2 and expand this circle to a circle with infi- 
nite radius. If the circle reaches the cut, it will circum- 
vent this cut along the real axis, resulting in two line 
integrals. In assuming that in the limit of an infinite 
radius the circle part of the integral vanishes, 
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This identity, known as dispersion relation, is valid not only for the used example but 
generally if H m (q 2 ) falls off sufficiently fast for \q 2 \ — > 0. However, this is not the case for 
the correlator function we have to deal with in our real application, 
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-DiscU(s), DiscU(s) = U(se- i7T ) - U(se i7T ) (8) 
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(it is convenient to write the correlator function in Euclidean space). In this case II(g 2 ) is 
singular. But we can redefine IT(g 2 ) by a subtracted quantity, 
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This subtraction method is known as momentum subtraction at a specified momentum 
square, in this case at the point q 2 = 0. Taking only the hadronic contribution to the 
correction, the spectral density is related to the relative hadronic cross section in e + e~ 
annihilations 

cr(e + e~ — > hadrons) 



R 
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by R(s) = 12tt 2 p(s). For low energies this relative cross section is determined by the pion 
form factor F n , 
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So the expression we deal with for the hadronic contribution to the fine structure correction 
is given by 
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The integrand consists of the spectral density and a factor which is called weight function. 



4 



Hadronic contribution to the anomalous magnetic moment 

The anomalous magnetic moment can be determined by 
looking at the first order QED correction to the muon- 
photon vertex with an self energy insertion for the loop 
photon. Using again the dispersion relation, we obtain 
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where the weight function is H(s) = K(s)/s with 
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Details about how to obtain this function are shown in Appendix A. The integral can be 
calculated and gives rise to 
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where 



1 
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In this case K(s) has a cut along the real axis so that the total singularity of the weight 
function H(s) = K(s)/s at the origin is stronger than in the case of the QED coupling. 
If we again only take the hadronic contribution, p(s) is meant to be the spectral density 
corresponding to the hadronic correlator function. 



Local and global duality 

The question may arise why we cannot exclude the experimental values from our consid- 
erations at all and take the relation in Eq. (|5p in order to calculate the spectral density 
from a correlator function calculated perturbatively. The reason is that there is actually 
an obstacle in using this relation. As depending on methods of functional analysis, the 
inverse of the dispersion relation is only valid if there are no poles encircled by the path 
in the complex plane, a condition which is necessary to obtain this relation. These poles 
can have their origin from weight functions in combination with the spectral density. This 
means that if there is such a weight function included in the integration of the spectral 
density, the inverse relation shown above is only valid globally and not locally. We call this 
global resp. local duality. Two paths out of this situation are shown in the following. 
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3 The polynomial adjustment method 



A singularity in the encirled part of the complex plane occurs for the two parameters we 
want to look at. Therefore we cannot fully keep out the experimental measurements from 
our considerations. There is, however, a way to include them in an optimal way. This is 
done by our method |12| which will be presented in the following section. In the meantime 



this method has also been used by other authors (see e.g. Ref. The method is based 

on the fact that we can use global duality when the weight function is non-singular. This 
is the case for a polynomial function. Therefore, we might mimic the weight function by a 
polynomial function obeying different conditions which I will explain later. By adding and 
subtracting this polynomial function Pn(s) of given order N to the weight function H(s), 
without any restrictions we obtain 



p(s)H(s)ds = f S " p{s) (H(s) - P N (s)) ds + p p{s)P N {s)ds (17) 

J S n J S„ 



where [s a , Sb] is any interval out of the total integration range. But because the second 
term has now a polynomial weight, we can use global duality to write 

H p(s)P N {s)ds = — r DiscU(s)P N (s)ds = 

Js a 2llZ Js a 

= u(-s)P N (s)ds + ^- I Tl(-s)P N (s)ds. (18) 

2m J\s\=sa 2m J\ s \=s b 

Thus this part can be represented by a difference of two circle integrals in the complex 
plane. On the other hand, the difference H(s) — Pn{s) suppresses the contribution of the 
first part. Our method consists thus of the following steps: 



• replacing p(s) in the first part of Eq. (|T7[) by the value of the experimentally measured 
total cross section R(s) (see e.g. Ref. |I||), 

• replacing the circle integral contribution to flavours at their threshold by zero, 

• in all other cases inserting the QCD perturbative and non-perturbative parts of II(— s) 
on the circle. 

These replacements can be seen as a concept within QCD sum rules. To obtain the best 
efficience of our method, we have to restrict the polynomial function by the following 
contraints: 

• The method of least squares should be used to mimic the weight 

• However, the degree N should not be higher than the order of the highest pertur- 
bative resp. non-perturbative contribution increased by one (this is a consequence 
of the Cauchy's theorem which is involved in the analytical integration of the circle 
integrals) . 
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• Especially for the low energy region, the polynomial function should vanish on the 
real axis to avoid instanton effects. 

• In regions where resonances occur, the polynomial function should fit the weight 
function to suppress those contributions which constitute the highest uncertainty of 
the experimental data. 

What are the ingredients for this mixture then? We have to look at what experimental 
data we can use and how the perturbative expansions for the correlator function are given. 



The experiment side 



Data from experiments are crucial in the low energy range and in threshold regions where 
perturbative QCD cannot be applied. Here combined data sets from various electron 
positron annihilation experiments are used M which are accomplished by recent BES mea- 



surements 14 



In addition the use of precise r-decay data from Ref. |L5| by isospin rotation promises 
to be a rewarding step in the low energy region. The vector r spectral functions are related 
to the isovector e + e~ cross sections for the corresponding hadronic states X by fO^ 



a I=1 (e + e 



Ana 2 



-Vj=l(T 



X-u T ). 



(19) 



vj=i(t — > X u T ) is obtained by dividing the normalized invariant mass-squared distribu- 
tion dNx~ /Nx-ds for a given hadronic mass a/s by the appropriate kinematic factor, 



Vj=l(T 



X~v T 
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(20) 



where \V u d\ = 0.9752 ± 0.0007 denotes the CKM weak mixing matrix element |16| and 
Sew = 1 + 5 ew — 1-0194 accounts for electroweak second order corrections [[HJ. The 
spectral functions are normalized by the ratio of the respective vector branching fraction 
B(t~ — > X~v T ) to the branching fraction of the electron channel B(r~ — > e~v e v T ) = 
17.79 ±0.04 ||. 



The theory side 



The two-point correlator |18| is given by 

i f(0\ 3 r(x)j^(0)\0)e^d*x 



-g a pq 2 + q a qp)~R(q 2 



(21) 



where we only included the isospin contribution / = 1, in contrast to corresponding con- 
siderations for the t decay. The scalar correlator function II(g 2 ) consists of perturbative 
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and non-perturbative contributions which we include to the extend we need them to keep 
the accuracy. For the perturbative contribution to the correlator we use a result given in 
Ref. [EDI. I only write down the first few terms, 
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with L = ln(/i 2 /g 2 ) while in Ref. [EJJ the expression is given up to the order 0(a 2 , m 12 /q 12 ). 
The number of active flavours is denoted by rif. For the zeroth order term in the m 2 /q 2 
expansion we have added higher order terms in a s , 
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with ko — 1, k\ — 1.63982 and ki = 6.37101. We have denoted the yet unknown con- 
stant term in the four- loop contribution by C3. Remark, however, that the constant non- 
logarithmic terms will not contribute to our calculations. The non-perturbative contribu- 
tions are given in ref. ETJ, 
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where we have included the /^-contribution arising from the unit operator. In this 
expression we used the SU(3) colour factors CV = 4/3, Ca = 3, and Tp = 1/2. For the 
coupling constant a s as well as for the running quark mass we use four-loop expression 
given in Refs. |22|, Ej|, pJJ even though in both cases the three-loop accuracy would already 
have been sufficient for the present application. 



The evaluation of the QED coupling 

The calculation for the running fine structure constant or QED coupling is published in 
Ref. [0. An important step in preparing for the application of our method is to select 
points s a and Sb for the limits of the integrals resp. the radii of the circles. Except for the 
threshold regions there seems to be a wide range for placing these points. Indeed we have 
shown that the method is fairly independent of this choice but is in some cases limited by 
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interval for y/s 


N 


data 
contr. 


contribution 

to AofiiCM*) 


error 
due to Aj^g 


[0.28 GeV, 3.1 GeV] 
[3.1 GeV, 9.46 GeV] 
[9.46 GeV, 30 GeV] 
[30 GeV, 40 GeV] 
y/s > 40 GeV 


1,2 
3,4 
3,4 
3,4 


24% 
0.3% 
1.1% 
0.15% 


(73.9 ± 1.1) x 10" 4 
(69.5 ±3.0) x 10" 4 
(71.6 ±0.5) x 10" 4 
(19.93 ±0.01) x 10" 4 
(42.67 ±0.09) x 10~ 4 


0.9 x 10" 4 
1.4 x 10" 4 
0.06 x 10" 4 
0.02 x 10" 4 


total range 






(277.6 ±3.2) x 10~ 4 


1.67 x 10" 4 



Table 1: Contributions of different energy intervals to oi\J /i (Mz)- Second column: choice 
of neighbouring pairs of the polynomial degree N. Third column: fraction of the contribu- 
tion of experimental data Fourth column: contribution to Aa^ d (Mz) with all errors 
included except for the systematic error due to the dependence on A^g which is separately 
listed in the fifth column. 

computational constraints as for instance the fact that matrices fail to be inverted in some 
special cases. For this reason it is not advisable to make the intervals too narrow. 

As an example for a threshold case I show the first contribution, 

f 1 R{s)H{s)ds = r R cxp {s){H(s) - P N (s))ds + 6m I U QCD {-s)P N {s)ds. (25) 

J so J sq J\s\=si 

Here we selected the range from the light flavour production threshold Sq = 4m 2 going to 
the next threshold marked by the mass of the ip, s\ = m\ ~ (3.1 GeV) 2 . In this case of 
course the inner circle integral did not contribute. The polynomial weight functions are 
compared to the weight function itself in Fig. |l|, the results for different polynomial degrees 
are shown in Fig. ^|. Further subdivisions can be found in Tab. |]. The part of the integral 
starting from S4 = (40 GeV) 2 up to infinity is done using local duality, i.e. by inserting the 
function R(s) obtained for perturbative QCD, 

(N c is the number of colours) into the second part of Eq. fllT|). 

To obtain the results shown in Tab. [I], we used the condensate values 
01 

{-LGG) = 0.04 ± 0.04 GeV 4 , a s (qq) 2 = (4 ± 4) x 10" 4 GeV 6 . (27) 

7T 

For the errors coming from the uncertainty of the QCD scale we take 

A^ = 380 ± 60 MeV (28) 
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The errors resulting from the uncertainty in the QCD scale in different energy intervals 
are clearly correlated and will have to be added linearly in the end. We also include the 
error of the strange quark mass in the light quark region which is taken as 

m s (lGeV) = 200 ±60 Me V (29) 

For the charm and bottom quark masses we use the values 

fh c (m c ) = 1.4 ±0.2 GeV, rh b (m b ) = 4.8 ± 0.3 GeV. (30) 

Summing up the contributions from the five flavours u, d, s, c, and b, our result for the 
hadronic contribution to the dispersion integral including the systematic error due to the 
dependence on Aj^g (column 5 in Table [I]) reads 

Aa^ d (M z ) = (277.6 ± 4.1) x 10" 4 . (31) 

In order to obtain the total result for a(Mz), we have to add the lepton and top contri- 
butions. Since we have nothing new to add to the calculation of these contributions we 
simply take the values from Ref. | |25| ], 

Aa\ iad (M z ) = (-0.70 ± 0.05) x 10" 4 , Aa icp (M z ) « 314.97 x 10" 4 . (32) 

Writing Aa(M z ) = Aa Xcp (M z ) + Aa had (M z ) our final result is 

a(Mz)- 1 = a(0) -1 (l - Aa(M z )) = 128.925 ± 0.056 (33) 

where we used a(O)" 1 = 137.036. 



The evaluation of the magnetic moment 

The procedure we choose for the anomalous magnetic moment of the muon differs from 
the one used for the QED coupling because of two reasons. As mentioned earlier, the 
singularity of the weight function at the origin on the one hand is much stronger and it 
therefore it would need more effort to approximate the weight by a polynomial in the lower 
energy region. On the other hand we can and should therefore make use of the fact that the 
r data between s = (0.28 GeV) 2 and si = (1.4 GeV) 2 were measured with a high precision, 
so they should influence the evaluation of the dispersion integral at full extent. The less 
precise e + e~ data from the above energy region, however, are worth being replaced by 
QCD expressions as much as possible. 

Our results are presented in Table 0. The dealing with error estimates is at some point 
different from the precious case of the QED coupling because of the different data sets we 
used. We only keep error estimates steming from the uncertainties for the quark masses 
and a systematic error over all ranges from the uncertainty of the parameter A^jg. In 
comparison to this, all other error estimates for the theory contribution of the "mixed" 
regions (like errors of vacuum expectation values) can be neglected. The value obtained is 
comparable with the predictions in Ref. |L3| and can therefore be seen as one of the most 
accurate on this field. It will be published hopefully soon [f26]]. 
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Figure 1: Weight function H(s) and polynomial approximations Pn{s) in the lowest energy 
interval 2m 7T < Vs < 3.1 GeV. 
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Figure 2: Comparison of the l.h.s. and r.h.s. of the sum rule given by Eq. fl25|) in the 
interval 0.28 GeV < Vs < 3.1 GeV. Dotted horizontal line: value of integrating the 
l.h.s. using experimental data including error bars [|TJ. The points give the values of the 
r.h.s. integration for various orders N of the polynomial approximation. Straight line 
interpolations between the points are for illustration only. The dashed lines indicate the 
error estimate of our calculation. 
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Fl 4 GeV 3 1 GeVl 
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(52 13 + 2 04) x 1 


-10 


nol vn nmia 1 mpt Vi nrl 


J ftp resonances 


(8.81 ±0.61) x 10" 


-10 


e + e~ annihilation data 


[3.1 GeV, 40 GeV] 


(22.13 ± 1.14) x 10 


-10 


e + e~ annihilation data 


T resonances 


(1.21 ±0.19) x 10" 


-10 


e + e~ annihilation data 


[40 GeV, oo] 


0.15 x 10- 


-10 


theory 


top quark contr. 


< 10 


-13 


theory 


whole range 


±1.83 x 10 


-10 


uncertainty from Aj^g 


hadronic contr. 


(694.1 ±7.0) x 10 


-10 





Table 2: The different contributions to the hadronic part of the anomalous magnetic 
moment a^ ad of the muon. 



4 Estimates for the next-to-leading order 

At this point we come back to the point whether the experimental measurements for the 
hadronic e + e~ channel really have to be taken into account. The question is to be taken 
serious because in this case it is at least strange to denote the estimate "theory estimate" . 
The method to use experimental data for a theory estimate is different from the (in this 
respect more preferable) method to use the experimental data in order to extract a few 
parameters like the locations of resonances and their widths. In order to see whether 
such an approach is possible, we use simple models to extract a single parameter from the 
given value of the anomalous magnetic moment of the muon. On the first sight this looks 
like a "zero-sum game", but using the model we then can estimate the next-to-leading 
order contribution to the anomalous magnetic moment which is found to be very close to 
data-based results in the literature. 



Integration in the Euclidean domain 

As mentioned earlier, the non-analytic structure of the complex plane does not allow for the 
calculation of the spectral density (given in Minkowskian domain) by using perturbation 
theory results (obtained in Euclidean domain). Therefore, at this point we return to the 
Euclidean domain. The transition from 

a 2f a \ 2 f°° K ( s ) / \ ; v f 1 z 2 (l-z)dz , nA , 

a u = An 2 - / — —p(s)ds, K(s) = / -= ) [ , = (34) 

\nj hm\ s HK ! ' V ; Jo z 2 + (l-z)s/m 2 v ; 

to the integral representation 

S = 4 ^ (^) 2 /o" W{t) (r nhad W) dt ( 35 ) 
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where 



W(t) 



Am 4 



y/t 2 + AmH (t + 2m 2 + y/t 2 + 4m 2 t) 
is calculated in Appendix B. Finally, an integration-by-parts leads to 

/ dU had (tY 



(36) 



An'l- 

,7T 



dt 



F(t)dt 



where 



F(t) 



W(t')dt' 



1 ft + 2m 2 - y/t 2 + AmH 



2m 4 



2 \t + 2m 2 + VWT&nHj (t + 2m 2 + Vt 2 + AmH 



2 ■ 



(37) 



(3* 



For high values of the Euclidean squared energy t the expression — rfil had (— t)/dt can be 
taken from perturbation theory while for the whole range different models can be used 
which are introduced in the following. However, we can take a short look on the weight 
function F(t). This function shown in Fig. |3]for m = m M = (769.9 ±0.8) MeV ]T6| supports 
to very high extend low energy contributions. This will become important in the following. 



,0.5 




t [GeV z 



Figure 3: The LO Euclidean weight function F(t) 



The step model 

The simplest model we can find for the e + e~ hadronic spectrum is the model of steps which 
appear when a new quark flavour is opened, i.e. when quark-antiquark pair production 
becomes kinematically possible at this specific energy. The model would simply read 

Pi(s) = N c Q 2 6(s - Ami) (39) 
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where N c = 3 is the number of colours and Q is the charge (in units of the elementary 
charge) of the quark flavour with mass tuq. Actually we need not separate between the 
three lightest flavours but instead can take them in combination, 

Pi ds = N C (Q 2 U + Q 2 d + Q 2 s )9(s - 4m e 2 ff ) = N c Ql s 6(s - 4m c 2 ff ), Q 2 ff = | (40) 

where m cS is the mass parameter which has to be adjusted in this model. For the step 
functions of the heavy quarks we use their masses from literature as cited before. Starting 
with this model spectral function, we obtain for the hadronic correlator function 

™ = 4 C a Kfr = - NAt Q *k = NA ln trfe) (41) 

and for it's derivative 

dUf s (t) _ N c Q 2 eS 



(42) 



dt t + Am 2 Q ' 
This actually makes sense, since perturbative QCD for high predicts 

_ dU^jt) = NcQtz f a s (-t) \ 

dt t \ it ) 1 1 

While for small values of t this expression remains finite. 
The threshold model 

One can of course use more sophisticated models. One of these is given by the fermionic 
spectral density 

P2(s) = N c Q 2 c ^l- —l[l + —iy (44) 
The corresponding hadronic correlator function is given by 

n 2 (t) = iV c Q c 2 ff {Q-3)^)-i} (45) 

where 

y^) = -^artanh (v ^)-l, z = J ^ (46) 



and 



dt iVc ^ff| t 2 + t 3 M t + 4m2 \\t + 



artanh 



4m 2 , 



(47) 



It figures out that the derivatives of II" ds (t) and II2 (t) coincide at the origin if we chose 
m eS /m q = y/5/2. 
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The resonance model 

Finally, one can also include the lowest lying resonance, namely the p meson resonance. 
A one-scale no-parameter model that satisfies the duality constraints from the operator 



product expansion is given by pTL 28, 29 



Pa(s) = N c Q 2 cS (2m 2 p 5(s - mj) + 9(s - 2m 2 p )) . (48) 
The derivative of the corresponding correlator function reads 



dt Weff \(t + m*) 2 ' t + 2ml 



NcQl s j—-^ + ■ (49) 



Alternatively, the narrow width p resonance can be replaced by a Breit-Wigner resonance, 
p 3 (s) = N c Q 2 cS (6(2m 2 p - s)p BW (s)6(s - Ami) + 0(* - 2mJ)) (50) 

where 

t \ 2m2 p ( T f m P \ (v\ 

PBw{) ~ ~ ll^FW+W ) ' 

the derivative of the correlator function is given by 

d Mt) _ N(1 2 ( f 2m ? PBw(s)ds 1 \ 

" ~ ! ^ cQcS [Li (s + ty + TT2^ P ) ■ { } 

The big surprize 

The spectral densities we have used are quite different. A glance on Fig. § convinces us 
that at least the spectral density of the third model is quite different from the two others. 
However, calculating the derivative of the correlator function for the three models, we 
obtain the result shown on the top of Fig. |5|. The two curves of model 1 and model 2 
coincide at the origin because we chose the relation m eS /m q = VE/2. We can even adjust 
m e ff/m q (and, therefore, m q ) in order to obtain a conicidence with the third curve at the 
origin. The surprizing result is that despite the fact that the spectral densities in Fig. |] 
look quite different, all curves for the correlator derivatives have nearly the same shape. 
But because the weight function F(t) supports contributions at low values of t this means 
that the method is independent of the choice for the modett The procedure, therefore, is as 
follows: we take the first model as the easiest one to obtain our results. The contributions 
to the anomalous magnetic moment of the muon for the heavy quark flavours can be 
calculated to be a^ ad (LO; heavy) = (71 ± 18) x 10~ n . The leading order result from the 
last chapter is reduced by this value to result in a^ ad (LO; light) = (6870 ±72) x 10 -11 . This 
value is used to adjust m e g. We obtain pi 



m eS = (201.7 ± 1.2) MeV. (53) 
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Figure 4: s-dependence of the spectral functions pi(s) of model 1 and P2(s) of model 2 
(upper diagram), as compared to the spectral function p-s(s) for model 3. We use m eS = 



201 MeV, m q = 179 MeV, and the central values m p = 769.9 MeV and T p = 150.2 MeV gB|. 
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Figure 5: The functions are the negative derivatives of the correlator function for the 
three different models. m cS = 200 MeV is used for the upper diagram and m eff = 205 MeV 
for the lower diagram. The parameter m q used for the model 2 is connected to m c g by 



m q = 2m eff / v / 5. The values m p = 769.9 MeV and T p = 150.2 MeV are taken from Ref. || 
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(a) (b) (c) (d) (e) 

Figure 6: LO and NLO corrections to the anomalous magnetic moment 

Next-to-leading order estimate 

As I mentioned earlier, up to this point the method seems to be a "zero-sum game". 
However, there is really makes sense if we want to calculate the next-to-leading order 
correction. The correction shown in Fig. 0(b) is the convolution of the spectral density 
with the two-loop kernel K^ 2 >(s), 

aJ ad (NLO(b)) = 4tt ( -J J ^ y } p{s)ds. (54) 



Using model 1, we obtain pTJ a^ ad (NLO(b)) = (-211 ± 5) x 10 _n . For the two bubbles 



in Fig. 0(c), the so-called double-bubble diagram, model 1 can again employed to give pT 
aJ^ ad (NLO(c)) = (106 ± 2) x 10~ n . Both results, as well as their sum are very close to the 



data-based results of the literature |32|. Finally, we add also the two next-to- leading order 



contributions, employing the four-point correlator, namely the light-by- light contribution 



(Fig. |^(d)) and the two photon Green function (Fig. 0(e)) (see the discussion in Ref. PIP ) 
to obtain our estimate 

aJ ad (NLO) = (85 ± 20) x HT 11 . (55) 
This estimate is still in good agreement with the difference 

a° xp (NLO) = (7165 ±151| cxp ±2.9| qed ±4|ew- (6941 ±70)) x 10- 11 = 

= (224 ± 189) x 10~ n (56) 

between the experimental value and the theory estimate for the LO correction. The im- 
provement of the accuracy of the experiment will show whether this can be kept or not. 



5 Conclusion 

I have presented results for the QED coupling a at the Z° pole and the anomalous magnetic 
moment a M of the muon, showing the main features of the polynomial adjustment method. 
In the second part I have used the LO theory estimate to adjust a mass parameter and 
calculate the NLO contributions. 
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How does the function K(s) appear? 

p'+ k 




(q = 0) 



If we look at the figure at the left hand side which 
is the diagram we have to consider, the correlator (in 
Minkowskian domain) shown as shaded bubble is given 
by 

IV(-A; 2 ) = (k p k a - k 2 g pa )U(-k 2 ). (57) 
For the scalar correlator function we take 



k 2 p(s)ds 



(58) 



The contribution for the diagram, therefore, is given by 
d D k 



x k 2 [9*p 



\ k a kp 



+ ft — m 
-ik 2 ) (V 



(-ief 4 ); 
k p k a \ 



+ ft — m 



k 2 



m-k 1 



—i 



d D k _ -f a (f + ft + m)^ (i + f + m)^ 
pU{p )- 



(27T) 



{{p> + k) 2 - m 2 ){{p + k) 2 - m 2 )k 2 



9a/3 



9a(3 



k a k/3 



X) 



kgkp 

k 2 



k 2 



Il(-k 2 ) 



(59) 



while the anomalous magnetic moment is the part of the diagram proportional to {p + p'Y 
for q 2 = (p — p') 2 = 0. The contraction with the second part of the effective transverse 
photon propagator simplifies the expression remarkably. Because of 



u(j/)ffl + ft + m) = u{p')[2p'k + k 2 + (-/ + m)lf) = u(p')(2p'k + k 2 ), 

{j> + ft + m)$ = (2pk + k 2 + + m)) u(p) = (2pk + k 2 )u(p) (60) 

(note that u{p'){j/>' — m) = {jf> — m)u{p) = 0), we obtain the two denominator factors 

ip' + k) 2 - m 2 = p' 2 + Ip'k + k 2 - m 2 = Ip'k + k 2 , 

(p + k) 2 -m 2 = p 2 + 2pk + k 2 -m 2 = 2pk + k 2 (61) 

(p 2 = p' 2 = m 2 ) which cancel out. However, the contribution figures out to be proportional 
to 7 M and therefore does not contribute to the anomalous magnetic moment. Therefore, 
we can contract with g a p and simplify the Dirac structure. 
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Simplification of the Dirac structure 

The simplification is done systematically using u{p' ){]/>' — m) = (jf> — m)u{p) = 0, 7 Q 7 Q = D 
and the usual relations for the Dirac matrices, 

= u(p')(2(p' + k) a + (-/ -lt + m) 1 a )Y(2(p + k) a + la (-^ -ft + m))u(p) = 

= u(p')(2(p' + kr-lt 1 a )Y(2(p + k) a - 1 J)u(p) = 

= A(pp' + (p + p')k + k 2 )u{p')^u{p) + u{p')^ a ^^u{p) + 

-2«g»V(/ + mu(p) - MpW + fth^(p) = 

= (4m 2 + 4(p + p')k + 4k 2 )u(p')^u(p) + (2 - D)u{p')$^$u{p) + 
-Ap ,fl u(p')ltu(p) + 2mu{p')^u{p) - 2k 2 u(p)^u(p) + 
-4,jfu(p')^u{p) + 2mu(p')^u(p) - 2k 2 u(p')^u(p) = ... 

= u(p) ((4m 2 + 4(p + + (£> - 2)A; 2 ) 7 M + 

-(4(p + p'f + 2 (£> - 2)fc")# + 4mF) M (j)) (62) 

where = g 2 = (p — p') 2 = 2m 2 — 2pp' =^> pp' = m 2 was used. Again, the first part is 
proportional to 7 M and does not contribute to the anomalous magnetic moment. 

Scalar, vector and tensor integrals 

Inserting the dispersion relation for Il(— k 2 ), the denominator factor k 2 is replaced by (s — 
k 2 ). The basic integral to deal with will be calculated by using Feynman parametrization, 

r d D k 1 

J (2it) D [k 2 + 2pk)(k 2 + 2p'k)(k 2 - s) 

d D k_ T ^ f 1 f 1 '* 1 dz x dz 2 



I 



(2tt) d v ' Jo Jo ( Zl (k 2 + 2pk) + z 2 (k 2 + 2p'k) + (1 - z x - z 2 )(k 2 - s)) 3 

f 1 A f 1Zl A f 1 

^ ' Jo Zl Jo 22 J {2n) D {k 2 + 2{z lP + z 2 p')k -(l-z x - z 2 )sf 
a ri~zi r dPk 1 
=: r(3 U dZl Jo dZ2 J (2^ (k 2 + 2Pk - M 2 f = 

rl rl—zi r d D fc 1 

= r(3)y o d Zl J o dz 2 J (27r)D — - — _ p2 _ M2)3 = 

= T(3)J q dz,J Q dz 2 J (2?r)D {k/2 _ p2 _ M2)3 = 
-iT(3-D/2) f 1 f 1 -^ 



( 47r )D/2 



/ dz x [ Zl dz 2 (P 2 + M 2 )^ 2 " 3 (63) 
Jo Jo 
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where for the last step a standard integral was used which is calculated by employing the 
Wick rotation. Note in addition that a shift k — > k' = k + P was performed. If we start 
with a vector integral including a factor k" in the numerator, one obtains 

r d D k' (k'-py 



7" = T(3) / dz! / dz2 / 
JO Jo J 



(2n) D (k' 2 -P 2 - M 2 ) 3 
<r(3-£>/2) _ /-i-*! 



[ d Zl [ 1 dz 2 P"(P 2 + M 2 ) D/2 ~ 3 , (64) 

JO JO 



the integrand proportional to k ,fl vanishes because the function is odd in k!" . For the tensor 
integral we obtain 

iv» = ~ lT j- 3 ~ D J 2 ^ f 1 dz x f 1 '* 1 dz 2 P»P v (P 2 + M 2 )^ 2 ' 3 + g^J =: I»» + g^J. (65) 
(47r) 1J / 2 Jo Jo 

Combining both ingredients 

If we combine the Dirac structure and the integrals, we obtain (up to general factors) 

- 4(p + p')^J x - 2(D - 2) lp I^ - 2(D - 2)PV + Ami". (66) 

Again, the term proportional to J does not contribute, we do not have to care about this 
integral in the following. Including all general factors, we obtain 



e'L' [6- U 2) f 

x iZ(p') (4(p + - 2(73 - 2)P"f - 4mP^) u(p) (P 2 + M 2 ) D/2 ~ 3 . (67) 
The integral is finite. Therefore, we can replace D = 4. Now we use 

u{p')fu{p) = u(p')(zi^ + z 2 f")u(p) = m(z 1 + z 2 )u(p')u(p) (68) 
and write M — (1 — z\ — z 2 )s and 

P» = z 1 p" + V = ^i + 22) (P + PT + ^1 - 22) (P "PT (69) 

where again only the first part counts for the anomalous magnetic moment. Finally, we 
can use z = z\ + z 2 as a new variable and obtain 

2 , ,w ,s „ , s Z" 1 , Z" 1 , 4m2 — 2mz 2 — 2mz 



A M = j^u(p')(p + pYu(p) J q d Zl J^ dz 



Izi m 2 z 2 + (1 — z)s 

e 2 r 1 . r* z(i-z) 



-u(p')(p + p'Yu{p) I dz I dz\ 



,2 



87r 2 m jo jo z 2 + (1 — z)s/m 

n u{p'){p + py u {p) f 1 9 z J~ z)dz , (70) 



2-kitl Jo z 2 + (1 — z)s/m? 

At this point the function K(s) appears. 
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Explicit form of K(s) 

The integral for K(s) can be calculated. First we can write 

f 1 z 2 (l-z)dz f 1 -z(z 2 + (1 - z)s/m 2 ) + (1 - z)zs/m 2 + z 2 

K(s) - 



z 2 + (l — z)s/m 2 Jo z 2 + (l — z)s/m 2 

f 1 , / s \ r 1 , r 1 zs/m 2 (2 — s/m 2 ) — s/m 2 (l — s/m 2 ) , 

... = - Z dz+ 1- — \ dz+ -L ^ ^— ^ f-A / ^. 

jo V m A ) Jo Jo z z + {l — z)s/m z 

The denominator in the integrand can be written as product 



z 2 + (l-z)^ = (z-z+)(z-z_) 



where 



s 2 / 4m 2 

z± = - — Al±v) = - , v — \l . 

2m 2 v ; Itv V s 



Therefore, we make the ansatz 

A + A_ (A+ + A_)z - (A + z_ + A_z + ) 
z — z + Z — Z- (z — z + )(z — Z-) 

and solve the system of equations 

A+ + A_ = A, (2 - A,) , A + z_ + A_z + = J-(\-J- 



m z \ m z / m A V m z 



to obtain finally 

(1 + v) 2 A (l-v) 2 

vA + = -j- J -, vA_ = ^ y 



Using 



we obtain 



l-z± = l 



{l-v) 2 ' ~ (1 + v) 2 ' 

2 l±v l-z± l±v 



1 =F v 1 =F v —z± 



K(s) = 1 ---±- 2+A+ ln( 1 -^±]+A_ln^- V - 



2 1 — v 2 \ —z + J \ —v. 

1 4 (1 + v) 2 ln ^ + ^ ■ (i-vf 



2 l-v 2 v(l-v) 2 V 2 J v(l + v) 

In using x = (1 — v)/(l + f ), we finally obtain 

w x 1 (1 + x) 2 1 + x / 1 \ 2 l + x_ / x 
K(s) = Tl -In + x 2 In' 



x x 2 (l — x) \l + xj l — x VI + x 

1 {1 + xf 1 + x ( 1 2 \ . „ , 2 1 + 2: 

+ — - a; 2 ln(l + x) + x 2 - lux 



2 x l — x Vx 2 / l — x 

which can be shown to be identical with the expression in the main text. 
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B Calculation details for the three models 

The main point here is to calculate the Euclidean analogon to the function K(s), i.e. W(t). 
We start with a reformulation of the denominator in the integrand of K(s), 

z 2 + (1 - z)— 2 = ( s + -^rrA (80) 

m z m l \ 1 — z J 

and conclude that 

t = m 2 z = (-t + Vt 2 + Am 2 t) (81) 

1 — z 2m 2 v ' 

is the substitution we have to use. We obtain 

if t + 2m 2 \ t + 2m 2 - Vt 2 + Am 2 t , 

dz = —1 H . at = , at, 

2m 2 V ^ 2 + 4m 2 t / 2m V* 2 + 4m 2 t 

z 2 = (t 2 + t 2 + AmH - 2tVt 2 + Am 2 t) = (t + 2m 2 - Vt 2 + Am 2 t) , 

Am 4 V / 2m 4 V / 

1-z = (t + 2m 2 - Vt 2 + Am 2 t) (82) 

2m z v ' 

and are lucky that a general factor occurs at all places. Moreover, for this general factor 
we can use 

(t + 2m 2 - Vt 2 + Am 2 t) (t + 2m 2 + Vt 2 + Am 2 t) = Am 4 (83) 

to switch freely from numerator to denominator and vice versa. For the denominator we 
obtain 

z 2 + (1 - z)^- = (t + 2m 2 - Vt 2 + Am 2 t) (s + t), (84) 
m 2 2m 4 v ' 

for the numerator we get 

2 t(t + 2m 2 -Vt 2 + Am 2 t) 3 ^ 

z (1 - z)dz = = dt, (85) 

8m 8 yt 2 + AmH 

therefore, finally 

Z" 1 2 2 (l-2)rf2 r 00 i(i + 2m 2 - Vt 2 + AmH) 2 

^ Jo z 2 + (1 - ^)s/m 2 ~ Jo 4m 4 v / t 2 T4^ 2 t(s + t) 



4m 4 t(it f°°tW(t)dt 



o Vt 2 + 4m 2 t(t + 2m 2 + Vt 2 + 4m 2 t) 2 (s + t) Jo s + t 



(86) 



Finally, we have 



a, = 4*>(°) 2 r ^lp(s)ds = An 2 («) 2 T f ^ P (s)dtds 



= Air 2 
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The Euclidean weight F(t) 

W(t) was not the final weight function we selected for our considerations. Instead, we 
performed an integration-by-parts, taking the derivative of H(t) and integrating W(t). 
The Euclidean weight function F(t) should be given by 

F(t) = J W(t')dt'. (88) 

We show this by checking that F(oo) = and F'{t) = —W(t) for the fiven 

Fit) = - ( t + ^-^ + ^t) = 2m 4 
U 2 \t + 2m 2 + Vt 2 + 4mH J (t + 2m 2 + Vt 2 + 4m 2 t) 2 ' 1 ' 

The first condition is granted. For the second one we take 



-4m 4 (l + (t + 2m 2 ) /y/t 2 + AmH) -Am 4 (t + 2m 2 + -Jt 2 + AmH) 

F (t) = - 



(t + 2m 2 + y/t 2 + 4m 2 t) 3 ^/W+lwH{t + 2m 2 + V* 2 + 4m 2 t) 3 

— zl 777 ^ 

= -w(t). (90) 



Vt 2 + 4m 2 t(t + 2m 2 + Vt 2 + 4m 2 t) 2 

The correlator for model 1 

The first calculation is quite easy and short, 

Ili(t) _ r°° pi(s)ds f°° 9(s-4m 2 s ) _ + f°° ds 



~ f / 



N cQls V s(s + t) 74ml S(S + f) hml s s{s + t) 

ds ds \ , / s 



r°° I ds ds \ / 

•A,,;, \ S S + t) 11 \ 



oft i + 4m 2 ff ' 



s + t 



4 ™cff 



(92) 



The correlator for model 2 

Starting with 

,r ^9 /! 4m 2 / 2m 2 \ 1 AT ^ 9 ,„ 9 . F 4m 2 4m 2 

p( s )=iV c g 2 ff yi-^h + -— ) =-N c Q 2 cS v(3-v 2 ), v = ^l-— : s = — - 

(93) 

we obtain 

n 2 (t) r°° ds F 4m 2 "/ 2m 2 \ r 1 v 2 (3-v 2 )dv 

nFQ^ ~ ~ J*™* s(s + *)V ~ ~V + ~ 4m 2 + t(l-t; 2 ) ~ 

_ /-i ^ 2 (3-^ 2 )rft; f 1 t; 2 (3-t; 2 )cfo 2 _ 4m 2 

~ 7o 4m 2 + t-^ 2 ~ 7o a 2 -v 2 ' + t ' 1 J 
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We continue with 

n 2 (t) f 1 v 2 (3-v 2 )dv y 1 2 2 y 1 2 2 r 1 dv 

= -l-^-V + ^-Vj!T^ = 4 + (« 2 -3)(a rta nh(I)-l). 

(95) 
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